Entanglement generation between unstable optically active qubits without photodetectors 
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We propose a robust deterministic sciieme to generate entanglement at liigh fidelity wittiout the need of 
photodetectors even for quantum bits, qubits, with extremely poor optically active states. Our protocol employs 
stimulated Raman adiabatic passage for population transfer without actually exciting the system. Furthermore, 
it is found to be effective even if the environmental decoherence rate is of the same order of magnitude as the 
atom-photon coupling frequency. Our scheme holds potential to solve entanglement generation problems, e.g., 
in distributed quantum computing. 
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One of the key challenges for practical quantum computing 
is scalability. Recently, an approach referred to as distributed 
quantum information processing has been suggested to solve 
this problem |[I}{6l. In this scheme, a scalable computer is 
constructed from a network of small devices, each compos- 
ing of a single or a few quantum bits, qubits. Importantly, the 
relatively long distance between the qubits renders it feasible 
to address the individual qubits and to suppress decoherence 
caused by unknown qubit-qubit interactions. To construct the 
network, inter-node entanglement is necessary, and many pro- 
posals of entanglement generation using a photon to mediate 
interactions between the qubits have been suggested [T'J?'.'?!. 
Due to weak interactions with its environment, a photon seems 
an ideal candidate for a flying qubit to generate such shared 
entanglement between computation nodes. 

However, such remote entanglement can be degraded by er- 
ror sources: imperfect photodetection and unstable optically 
excited states. Although there are many proposals to perform 
entanglement generation, in most cases, the use of photode- 
tectors for measuring the emitted photons is inevitable, and 
imperfections of detectors causes a significant loss of the fi- 
delity. The first experimental realization to perform entangle- 
ment generation between macroscopically distant atoms has 
been reported in Ref. fS]. The fidelity in this probabilistic 
method is only 0.63 after post selection, limited mainly by the 
dark counts of the photodetectors. In general, it is usually nec- 
essary to involve optical transitions into an excited state of the 
qubit for the entanglement generation. Such an excited state 
is prone to decoherence due to the strong coupling with the 
environment ll9llT0l. 

A protocol, in which the fidelity of the entanglement gen- 
eration does not depend on the imperfections of the photode- 
tectors, has been presented in Ref. |7 1. However, this scheme 
requires optical transitions and hence becomes sensitive to the 
decoherence of the optically excited state. Although other 
protocols have been suggested recently in order to overcome 
this type of decoherence ifTTl [121 . the protocols still require 
photodetectors to generate the remote entanglement. Conse- 
quently, they are vulnerable to the imperfect photon detec- 
tions. In this Letter, we introduce a protocol, in which the 
eventual fidelity is not hindered by either of these effects. Sur- 
prisingly, without the use of photodetectors, this protocol still 
achieves a high fidelity entanglement even when the environ- 



mental decoherence rate at the excited state is as large as the 
atom-photon coupling frequency. 

The basic idea of our scheme is to utilize the concept 
of single-particle entanglement suggested by van Enk fl3|. 
Suppose that one uses a half-plated mirror to split a single 
photon into two paths, and on each path there is an atom pre- 
pared in its ground state. The photon on each path is fo- 
cused onto the atom to be absorbed. For simplicity, let us 
make the assumption that the absorption probability is unity. 
Although one of the atoms will be excited, it is not possi- 
ble to know which atom is excited and therefore a Bell state 
l*i^^> = ^(|G)L|^)fl + |£;)L|G)fl) can be generated d 
where \G)i and \E)i {i = L, R) denote the ground state and 
the excited state of the i th atom {i ~ L, R), respectively. 

However, there are several potential sources of errors to 
challenge this simple scheme. Firstly, the excited state is 
usually prone to decoherence lfT0l[T4ll . Secondly, a low ab- 
sorption probability is a serious source of error in this sim- 
ple protocol, since the interaction between the photon and 
atom in a free space is weak. To overcome these difficulties, 
our scheme involves a qubit in a cavity forming a lambda- 
system with a ground state \G), an exited state \E), and a 
metastable state \M) as shown in Fig. [T] The state \M) 
is coupled to the excited state \E) through the cavity cou- 
pling strength g and we induce Rabi oscillations between 
the states \G) and \E) by a laser Initially, the state is pre- 
pared to be \M)L\M)ii\\ac) where |vac) denotes a vacuum 
state for the cavity photons. A single photon split by a half 
mirror is focused onto the two cavities which lie symmetri- 
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cally. As a result, we have |Af)i|A/)/j4=(a}^ -r- uj^j 

where a\ {i — L, R) denotes a creation operator of a cav- 
ity mode at the ith atom. Importantly, by ramping adiabat- 
ically off the classical field at each cavity, the population of 
the state |M)a^^|vac) can be transferred to the state |G)|vac) 
essentially without populating the excited state while the state 
|Af)|vac) remains intact lil5J . Thus the state evolves into a 
state -^{\G) l\M) R + |Af)i|G)fl)|vac). Similar adiabatic 
transfer between a photon excitation and a single atom in a 
cavity has been demonstrated experimentally in Ref fT6\. We 
stress that, in our scheme, neither a use of photodetectors nor 
an optical transition to the excited state is necessary, and there- 
fore this protocol is robust against typical errors caused by the 



imperfections. 
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FIG. 1 : Schematic picture of an apparatus for the entanglement gen- 
eration. A half-plated mirror splits a single photon into two paths, 
and at the end of each path a qubit is confined in a cavity. Each 
atom has a A-type energy level structure, i.e., the ground state |G), 
metastable state \M), and excited state \E). The state \M) is cou- 
pled to the excited state \E) with the cavity coupling strength g, and 
the state \G) is coupled to the excited state j_B) by a classical driving 
field A. 

In the above description, we assumed that the ramp of the 
classical driving field is so slow that an optical transition to the 
excited state does not occur. Furthermore, the photon loss dur- 
ing this process is negligible with a high-Q cavity. However, 
there is a non-zero probability to excite the state for a finite 
operation time, which reduces the fidelity of the entanglement 
generation. In order to include this effect, it is necessary to 
present a more detailed analysis. 

Let us consider first the case of a single-qubit system to 
study the effect of the decoherence of the excited state and 
the non-adiabaticity. Below, we generalize to the actual two- 
qubit case for entanglement generation. The Hamiltonian of 
the system is given by 



Hsyfi{t) — Hs + Hi, 



Hs 
Hi 



= e|e)(e| + uJcO, a, 



g(|l)(e| + |e)(l|) + A,cosM(|e)(0| + |0)(e 



where we have defined |0) = |G)|vac), |1) = |M)a^'|vac), 
|e) — |i?)|vac), and set h = 1 for simplicity. Here, At is 
the amplitude of the time dependent driving field, uj is its fre- 
quency, uic is the frequency of the cavity mode, e denotes the 
energy of the atomic transition, and g is the coupling constant 
of the standard Jaynes-Cummings model. For simplicity, we 
assume that the driving field, cavity mode, and atomic transi- 
tion are resonant so that to = ujc = £■ Since the classical field 
is ramped adiabatically off, the time dependent strength of the 
field can be represented as Aj = A cos^ (at) for < t < ^ 
where - denotes a time scale of the variation of the Hamil- 

a 

tonian. In a rotating frame with the laser frequency, which is 



characterized by a unitary operation e"^*!*^^^"! fTTJ, we obtain 

Hs 0, 6(|0)(0| + |e)(e| + |l)(l|), 

Hi c. 5(|l)(e| + |e)(l|) + ^^°^(|e)(0| + |0)(e|), 

where we have employed the rotating wave approximation. A 
convenient basis to describe the dynamics of this system is the 
adiabatic basis, defined as Hsys{t)\En{t)) — En{t)\En{t)). 
Also, in an appropriate frame to rescale the energy, we can 
assume e = without loss of generality. 

Since we have only one excitation in the system, the eigen- 
states become 



\Eo) 
\E,) 
\E,) 



Va 



cos'»(at)+4g 
i=Acos^(at) 



%/2 



y/4g^+A^ cos* {at) 



{-2g\0)+Acos^at)\l)) (1) 



v/2 



'(at) 



^4g2+A2cos<i(ai) 



A cos^ (at) 



^4g2+A2cos'i(at) 



The eigenvalues are £^{1) = 0, £'g(t) 

and E~,{t) = V^9'+^'^-o-Hat) ^ ^^ ^^ ^^j^ ^^ ^^^^^ ^^^ 

joint effect of decoherence and the adiabatic evolution, we 
aim to map the dynamical system into a time-independent 
one llT8l[T9l . To this end, we define a transformed system den- 
sity operator p^ys = D^ PsysD, the time evolution of which is 
governed by the effective Hamiltonian iJgff = H{t)+w. Here 
the Hamiltonian iJ(t) = D^H{t)D == Efc=o.e,i ^fcl'^)(^l 
is diagonal in the time-independent basis {|fc)}, the op- 



erator D = Efc=o,ea l^fc)(fc| 

2'\/2ig-a-A cos(ai) sin(ai) 



and w 



= -^D^'i 



4,^ +A^W(a.) (I0)(l| + I0>(e| - Il)(0| - Ie>(0|). 

The remaining time dependence is manifested in the correc- 
tion term w which tends to excite the system. There are a 
few possible strategies to treat this correction term. The sim- 
plest scheme is to disregard the correction term completely, 
which can be valid only in the adiabatic limit. To include the 
non-adiabatic correction to the lowest order, it is possible to 
perform another transformation coiTesponding to the super- 
adiabatic basis |20| as we will do in this paper. Namely, we 
diagonalize the Hamiltonian H{t) + w using the first-order 
perturbation theory on the correction term 1201 . Note that 
there is no energy shift in this order of the perturbation theory. 
The approximate eigenstates of the effective Hamiltonian in 
the untransformed system, i.e, eigenstates of D{H{t)+w)D^ , 
are expressed with the help of Eqs. ([T])-([3]) as 
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e)(6) 



4y2igaA cos{at) sin(ai) A cos (at) 



yt = 



and 



zt = 



[4g2+A2 cos" (at)] 5 """ ^4g2+A2 cos''(at) 

'^ As long as the adiabatic condition 



^4g2+A2cos4(at) 



|Bg_Eil ^ |Bg_Ee| ^ J jg satisfied, an initial state |£'g(0)) re- 
mains in the state \E={t)) during the adiabatic process ||2TI . 
Since we have |-B=(0)) « ll) for f < 1 and \E~J^)) = 
— 10), the population transfers from the state |1) to the state 
|0) when the driving field is ramped off adiabatic ally. In 
Eq. (|4]), ^/2xt denotes the amplitude of the excited state yield- 
ing P - 2ItJ2 ~ 6ig^ a^ A^cos^ (at) sin^ (at) . . „„„;.„ J. 

ing -Te — z|Xt| _ [4g2+A2cosi(at)]3 lor me exciieu 
state probability. Since we assume that g ^ A, the maximum 
probability is obtained at a point where sin(at) w 1. With 
this approximation, the maximum excited-state probability is 

given by ^j^^-r for t = ^ arccos \ -jtz- ^^ conclude 
that the slow ramp rate of the driving field prevents the adia- 
batically evolving state to have projection on the excited state, 
hence protecting the system from the decoherence of the ex- 
cited state. 

However, we did not take the decoherence processes explic- 
itly into account above. Especially, the optically excited state 
|e) is usually coupled strongly with the environment causing 
decoherence of the quantum states |9, 10 1. Although the pop- 
ulation of the excited state should be small in our scheme be- 
cause of the slow variation, we study how the noise degrades 
the coherence. Therefore, we need to derive a Markovian mas- 
ter equation. The total Hamiltonian of the system and the 
environment is represented as H ~ -ffsys(^) + ^im + ^env 
where iJenv denotes a bath operator and Hint denotes the 
coupling between the system and the environment. Since 
we consider the noise at the excited state, the coupling be- 
tween the system and the environment can be represented 
as iJint = X\e){e\ (E) X where X is the environmental op- 
erator. Equations (|4]l-<|6]l provide us conveniently the trans- 
formation operator corresponding to the superadiabatic ba- 
sis Ds ~ J2k=o e i\-^l)i^\ which yields the transformed 
density operator p^y^ = Dlpf^y^Dg. In this approximation, 
the time derivative of Dg is neglected and the effective to- 




FIG. 2: (Color Online) Population of the state |0) (green ascend- 
ing curve), |e) (red almost constant curve), and 1 1) (blue descending 
curve) during the adiabatic transfer under the effect of decoherence 
at the excited state |e). We have employed parameters a — 0.01 x g, 
A — 50 X g, and 7 = 0.1 x g. The results are independent of the 
value of g. The population of state 1 1) is transferred to the target state 
|0) with fidelity 0.992 when the external field is ramped off. 



To estimate the effect of dephasing of the adiabatically 
evolving state due to the decoherence at the excited state, 
we consider an untransformed initial state jV'o) — ^l-E-fi) + 
4=|r) where \r) denotes a reference state which is not cou- 
pled with any other states. Here, we assume that, due to the 
strong driving field ^ ^ 1, the effect of an imperfect initial- 
ization is negligible. Since we are interested in an early time 
stage of the decoherence process, we approximate 



kysj{t)^Dl\^o){'^o\Ds 



dt'[Li{t'),[Li{t'),~p,y,j{0)]], 



where we substitute Psya.iit) with /5sysj(0) in the inte- 
grand. We have employed the interaction picture defined for 
tal Hamiltonian becomes H = i?sy.s + ^int + ^env, where the operators as 0/(i) = Uld{t)Us{t), where Us{t) = 



V2xt\Q) 



E,=o.i e^^b')01. ^int = 10) (01 ® X, and 
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1) + 0(|xtp). Under the assumption 



of a white noise spectrum of the environment, we integrate the 
von Neumann equation for the total density matrix, trace out 
the environment, and arrive at the master equation 



e 'Jo "** ^sys(t ) jg (jjg system time-evolution operator. 
Thus, the fidelity can be calculated as follows: 

F = sm)\p.y^,s{t)m))s 



dps. 



dt 



i[Hsys,PsyM] - -^[L, [L,p,y,] 



(7) 



where 7 denotes the decoherence rate and L denotes a Lind- 

blad operator defined as L = 10) (01- Note that this noise 
operator can cause unwanted transitions from the adiabati- 
cally evolving state \E~){E~^\ « psys to the other states. We 
solve this master equation numerically and show the popu- 
lation of each state Pj{t) = {j\psys{t)\j) ij = 0,1, e) in 
Fig. I2] Starting from the state /Osys(O) — |1)(1|, we have 
achieved Po{j^) = 0.992 for a = 0.01 x g, A = 50 x 5, 
and 7 = 0.1 X g, which shows almost perfect transfer from 
the initial state |1)(1| to the target state |0)(0| even under the 
effect of decoherence at the excited state. 



(8) 



where 1-0(^)5 denotes the transformed adiabatically evolving 
state at a time t for vanishing decoherence. Since we have 
^ <C 1, 1 2;* I is negligible except for a time region satisfy- 
ing < A cos^ (at) < eg, where c is some constant. In this 
time region, we have < 7? — t < -\/^, where n is a 

^ ' 2a a V A' 

constant that depends only on c, not on a. Hence we obtain 

1-F W 7/f__ r^lXfA^dt' = ^7 f!^ „ r^P.dt' < 



\xt>rdt' = i7/r_ 

2a a V A 2a 

f^7^/f s^. Therefore we obtain 1-F = 0{a). As ex- 
pected, the harmful effect of the decoherence decreases with 
a. 

Let us generalize the above results to the two-qubit 
case. Since we consider independent decoherence processes 



for the qubits, the master equation for the two-qubit sys- 



tem IS 



dt 



- E 



j=L,B 



n-ffsys,Psys] — 2 [-^i i [^J 7 Psys]]- 

We solve the master equation numerically using the ini- 
tial state \M) l\M) f(-y={a]^ + a)j.)|vac), and show the fi- 

dehty (■(/'Beii|/0sys(5^)|V'Beii) in Fig. [5] Here, |i/'Beii) = 
-^{\M)l\G)r + |G')L|A^)fl)|vac). For reference, we show 
also the fidelity with neither rotating wave approximation nor 
transformation corresponding to the superadiabatic basis. The 
two methods give similar results which verifies that our en- 
tanglement generation takes place at high fidelity independent 
of the theoretical approach. Surprisingly, even when the envi- 
ronmental decoherence rate 7 is as large as the atom-photon 
coupling frequency g, the protocol still generates high-fidelity 
entanglement. 
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FIG. 3: Fidelity of the proposed entanglement generation scheme as 
a function of the decoherence rate of the environment, 7, with (top 
curve) and without (bottom curve) rotating wave approximation. We 
have fixed the parameters asa^O.Olxg, and A — 70 x g. 

Photon loss is one of the major problems in the most of 
previous entanglement generation protocols H] [3] ID 12214241 . 
In our case, it can also degrade the fidelity since the pho- 
ton might be lost before coming in the cavity, while in or 
entering the cavity, or the pulse generated by the single- 
photon source may fail to provide a photon 1251 . The 
state with photon loss can be described by plr = (1 — 
Pioss)|V'Beii)(V'Beii| + Pioss\M) l\M) R l{M\ m{M\ even for a 
perfect adiabatic transfer where Pios, is a probability to lose 
the photon. Fortunately, this is of the same form as the re- 
source state considered in Refs Q |26l l27l . If we have an 
ancillary qubit near the actual qubit at each location, we can 
perform an efficient two round distillation protocol as follows: 
It is possible to utilize the state p^R as a resource to perform 
the parity projection on the ancillary qubits. Although this 
parity projection is imperfect due to the photon loss, by gen- 
erating the state plr again through another adiabatic trans- 
fer, one can utilize the second state for performing the second 
parity projection on the ancillary qubits in order to obtain a 
perfect Bell state. The measurement results of the two par- 
ity projections let us know whether entanglement between the 
ancillary qubits is generated or not, and the success proba- 



bility is calculated P^ = ^^ ^'°"-' ll26l . Note that the parity 
projection is one of the most commonly proposed entangle- 
ment generation methods |l2l|7]|26| to make a two dimensional 
cluster state for quantum computation |28|. Moreover, the 
joint effect of the photon loss and decoherence in this distil- 
lation protocol has been studied |27|, and has been shown to 
be negligible for weak decoherence. These results show that 
this distillation protocol makes our scheme robust against the 
photon loss by using ancillary qubits near the optically active 
qubits. Such ancillary qubits are available with some species 
of nanostructures such as the nitrogen-vacancy centers in di- 
amond I 



In conclusion, we propose a deterministic scheme to gen- 
erate entanglement between unstable optically active qubits. 
Our method is designed to function without the need for pho- 
todetectors which are typically major sources of error. With 
the help of adiabatic transfer, we entangle the qubits with- 
out exciting the unstable states, which renders our proposal 
extremely robust against decoherence. The authors thank T. 
Close and B. Lovett for useful discussions. We acknowledge 
Academy of Finland, Emil Aaltonen Foundation, and Centre 
for International Mobility for financial support. 
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